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1 INTRODUCTION 

Weak currents can be classified as either first- or second-class depending on the J^^ of the decay 
current H], where G-parity is a combination of charge conjugation and an isospin rotation, G = 

Ce^^^, and is a multiplicative quantum number. In the Standard Model, first-class currents (FCC), 
where PG{—iy = +1 (J^*^ = 0++, , l"*"", 1^+, . . .), are expected to dominate decays while 
second-class currents (SCC), where PG(— 1)"' = — 1 (J^"^ = 0+^,0^+, 1++, 1 , . . .)/ are expected 
to be small and to vanish in the limit of perfect isospin symmetry. An example of such a decay 
is T~ -^ uj-n~v^ which is expected to proceed through FCC mediated by the p resonance. This 
decay may also potentially proceed through SCC, such as 6i(1235) |2J with r^ -^ h^Vr ^ ujtt~Vt, 
producing final state particles with J^*^ = 1++ and ^ . 

Since the decay h^ -^ ujir" occurs through S- and D-waves, as compared to a P-wave for FCC, 
different polarizations of u spin result in different angular distributions of the final state particles. 
The expected distributions of cos 9^t^ for all possible spin-parity states of the final state particles are 
listed in Table [H where 6^.,^ is the angle between the normal to the co decay plane and the direction 
of the remaining vr in the lo rest frame, as shown in Figure [H The existing measurement of the 
angular distribution of r^ —>■ ujTT~Vr is consistent with having only P-wave contribution, and the 

present limit is 5.4% for the ratio of SCC to FCC contributions, iV(';;Sn-vector current) / ^(vector current)' 

at 90% confidence level ^. This paper presents a search for SCC in r^ -^ ojii^Vr decays with 
uj — > vr+Tr^vr'^ by studying the angular distributions of final state particles. 

n 

A 




Figure 1: Illustration of the angle Oi^.,^: the angle between the normal to the lo decay plane and the 
direction of the remaining vr in the ui rest frame. 



2 THE BaBar detector AND DATASET 



This analysis is based on data recorded by the BaBm. detector [4] at the PEP-II asymmetric-energy 
e+e" storage rings operated at the Stanford Linear Accelerator Center. The data sample consists 
of 347.3 fb^ recorded at the center-of-mass energy of 10.58 GeV. With a cross section for r pairs 
of Urr = (0.919 ± 0.003) nb [5J, this data sample contains nearly 320 million pairs of tau decays. 

The EaEar detector is described in detail in Ref. [4J. Charged-particle momenta are measured 
with a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) inside 
a 1.5-T superconducting solenoidal magnet. An electromagnetic calorimeter (EMC) consisting of 



''charge-conjugate reactions are implied throughout this paper. 



Table 1: Expected angular distributions, Fl{cos 9^.,^), for possible spin-parity states in the decay of 
T" -^ uTT^Vr- L is the orbital angular momentum. 



J^ 


L 


ptVC 


^^■^■^'(cos^^.) 


1- 
1+ 
1+ 
0- 


1 

2 
1 


F^vu ^ (1 _ cos^e^^) 


F^c^ oc 1 

Fi^"^ oc (1 + 3 cos^e^,) 

Ff ^^ oc cos^e^. 



6580 CsI(Tl) crystals is used to identify electrons and photons. A ring-imaging Cherenkov detector 
is used to identify charged hadrons, in combination with ionization energy loss measurements 
(dE/dx) in the SVT and the DCH. Muons are identified by an instrumented magnetic-flux return 
(IFR). 

Monte Carlo (MC) simulations are used to estimate the signal efficiencies and background con- 
tamination. The production of r pairs is simulated with the KK2 f generator [6J, and the decays of 
the r lepton are modeled w^ith Tauola [7]. Continuum qq events are simulated using JETSET IH. 
Final state radiative effects are simulated for all decays using Photo s HI. The detector response is 
simulated with GEANT 4 ||T0|, and the simulated events are then reconstructed in the same manner 
as data. 



3 ANALYSIS METHOD 



Since r pairs are produced back-to-back in the e+e" center-of-mass frame, each event is divided 
into two hemispheres according to the thrust axis [TT], calculated using all reconstructed charged 
particle tracks. Candidate events in this analysis are required to have a "1-3 topology," where one 
track is in one hemisphere (tag hemisphere) and three tracks are in the other hemisphere (signal 
hemisphere). Events with four well-reconstructed tracks and zero net charge are selected. The 
polar angles, in the laboratory frame, of all four tracks and the neutrals used in vr'^ reconstruction 
are required to be within the calorimeter acceptance range. Events are rejected if the invariant 
mass of pairs of oppositely charged tracks, assuming electron mass hypotheses, is less than 90 
MeV/c^, as these tracks are likely to be from photon conversions in the detector material. 

The charged particle found in the tag hemisphere must be either an electron or a muon can- 
didate. Electrons are identified using the ratio of calorimeter energy to track momentum (E/p), 
the shape of the shower in the calorimeter, and dE/dx. Muons are identified by hits in the IFR 
and small energy deposits in the calorimeter consistent with expectation for a minimum-ionizing 
particle. Muons with momentum less than 0.5 GeV/c cannot be identified in this manner as they 
do not penetrate far enough into the IFR. Charged particles found in the signal hemisphere must 
be identified as pion candidates using dE/dx. The tt^ candidates are reconstructed from two sep- 
arate EMC clusters with energies above 100 MeV that are not not associated with charged tracks 
and are required to have invariant masses between 100 and 160 MeV/c^. Events are required to 
have a single vr'' in the signal hemisphere. The r candidates are reconstructed in the signal hemi- 
sphere using the three tracks and the n^ candidate, and the invariant mass of the r candidate, 
?Ti(27r~7r+7r'^) is required to be less than the nominal mass of the r lepton, rrir = 1.777 GeV/c^ lfT2| . 
After the event selection process, from the MC it is found that 14% of the events remaining are 
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Figure 2: a; candidate mass spectra for selected events in data and expected Monte Carlo back- 
ground (colored histograms). The background histograms do not include the non-resonant 
T~ -^ 2tt^ TT^ n^ Vj- decays. The signal (S) and sideband (SB) regions are indicated in the figure. 



r-pair events that do not contain a r 

For each selected event with 7n(27r^7r^ 



2tt TT^TT^Ur decay, and 1.3% are e~^e 



■ events. 



vr") < nir two CO candidates are reconstructed from 



vr+vr^vr combinations. The mass of the u> candidates, rn,(7r+7r^7r ), is required to be between 670 
MeV/c^ and 890 MeV/c^; within this range, the signal region is defined between 760 MeV/c^ and 
800 MeV/c^ with mass regions of width 60 MeV/c^ on each side of the peak used as sideband 
regions for background studies, as shown in Figure |2l For each m candidate in the signal region, 
the angle Oi^j^ is calculated and is used in the SCC measurement, after background subtraction. 

There are three background types to be considered in this analysis. The first type is com- 
binatoric background, which is expected to have an angular distribution that is independent of 
m{'iT^ TT^ TT^) , and is thus subtracted from the signal region using the sideband regions. The num- 
ber of combinatoric events lying w^ithin the signal region is obtained by fitting the m(7r''"7r~7r'^) 
spectrum with a smeared relativistic Breit-Wigner for the lo resonance and a polynomial for the 
combinatoric background. The polynomial is integrated over the signal region to find the number 
of continuum events in the signal region. The second type of background comes from e^e~ — > qq 
events that contain lo —>■ vr+Tr^vr'^ decays. While the event selection process significantly reduces 
the number of qq events, approximately 0.3% of the events in the signal region are expected to 
be of qq origin. This type of background is studied using events with m{27r^TT^TT^) well above 
the r mass (> 2.1 GeV/c^). In this region, where all events are considered to be of qq origin, a 
comparison of the numbers of events in MC and data is used to obtain a scaling factor for the qq 
background events. 

After subtracting combinatoric and qq background events, approximately 4.6% of the remain- 
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ing uj candidates in the signal region are expected to be background events from non-signal r 
decays. The dominant of these, comprising 99% of these background events, is t~ —f lutt^tt^Ut-, 
where one vr'^ has not been reconstructed. The decay t" — > Lun^ir^Ur has not been well measured 
and is incorrectly modeled in the MC. Both the branching fraction and decay angular distribution 
need to be corrected, as shown in Figures [Sfa) and |3tb). To correct for the differences between 
data and MC, events with an additional vr*^ candidate in the signal hemisphere are selected, using 
the same cuts discussed above. Using these events, the MC branching fraction of r^ —>■ utt^tt^Vt 
is corrected by comparing the numbers of fitted lo candidates in data and MC. The fit function 
used for this is a smeared relativistic Breit-Wigner with a polynomial background. The branching 
fraction obtained using this correction technique is found to be consistent with existing measure- 
ments ||T2| . To correct the angular distribution of t~ —>■ ujiT~'K^Vr in the signal region, backgrounds, 
consisting of combinatorics, qq events and r~ — > ujir^Vr decays, are subtracted from the two vr*^ 
data sample, and the remaining cos 0t^,r distribution, shown in Figure lUc), is used to correct the 
r^ -^ unr^Tr^UT distribution in the MC. 
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Figure 3: (a) 7n{n^ tt^ tt^) and (b) cos^^^^ distributions when requiring an additional vr'^ in the 
signal hemisphere, before background subtraction. These are used to correct the r^ 
MC. (c) The cos 6*^^^ distribution obtained from data after subtracting background. 
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To account for any variation in efficiency as a function of cos 9^^^, the generated and recon- 
structed MC cos0(^7r distributions of r^ -^ Loir^Ur in the signal region are compared. The ratio 
of the two distributions, shown in Figure IH is used as an efficiency function to correct the back- 
ground subtracted data. 
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Figure 4: Efficiency as a function of cos 9^^-,^ obtained from r ^ 27r 7r~^7:^VT- MC. 



4 PHYSICS RESULTS 



After subtracting background events and applying efficiency corrections, a binned fit to the re- 
maining cos 9i^Ty distribution is carried out using 



FCC/ 



Ficos9^^) =Nx [eFo*^^(cos0^O + (1 - e)F^^^(cose^^)] 



(1) 



where A^ is a normalization factor, the parameter e is the fraction of r — > wvr Ur decays that 
proceed through SCC, and F^^'~^ and F^^^ are normalized angular functions described in Ta- 
ble [D The parameter e is related to A^(non-vector current) / ^(vector current) by the equation e/(l - e) = 
^(nSn-vector current) /^(vector current)- I^ EqE Only F^^^ is used for the function describing the SCC 
contribution since the shape of this function gives the most conservative (largest) estimate of e. 

This method is tested by adding various amounts of ^-wave decays to the standard MC and 
fitting for the levels of SCC. As the results of this test indicate, as shown in Table |2l in all cases 
the measured fractions of SCC, e, are consistent with the fractions added to the MC. The errors 
listed in Table|2]are not necessarily indicative of the expected uncertainties in the data; they do not 
contain systematic uncertainties, and statistical correlations exist among the MC samples used in 
the studies. 

The largest contributions to systematic uncertainties on e are scaling and modeling of the MC 
background. The correction applied to the branching fraction of r^ -^ uj-k^tt^Vt has an error 
associated with it, determined by the available statistics. The correction factor is adjusted by ±lcr 
to obtain the uncertainty in e while the errors associated with correcting the angular distribution 
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Table 2: Test of SCC measurements in MC with various amounts of 5-wave decays added. 



Fraction of L = r ^^ urn Ur decays added 


e 


none 
1% 

2% 


(0.11 ±0.17)% 
(1.10 ±0.17)% 
(2.09 ±0.17)% 



are folded into the statistical uncertainty. In addition, there are r decays that may be present in the 
final event sample but which are not simulated in the MC. The largest of these are expected to be 
r~ -^ uK'^Vr, T^ — > LUTT^lir^Vr and r^ —f uj2tt~ 7r~^ Ur decays, which when combined can add up to 
0.2% of the final event sample. Since the effect that these decays have on the angular distribution 
is unknown, the extreme cases are taken to obtain the uncertainty. These cases correspond to 
these decays having either entirely 1 — cos^ Ounr or entirely cos^ Oujn distributions. The scaling of 
qq events can also affect the measurement of e, and the uncertainty is obtained by adjusting the 
scaling factor by ±lcr. These systematic uncertainties are summarized in Tabled 

Table 3: Summary of systematic uncertainties on e 



Source 


Uncertainty (a^) 


B (r^ —>■ ujtt~tt^Ut-) 

un-simulated r decays 

qq scaling 


±0.0007 

+0.0000 
-0.0055 

±0.0001 


Total 


+0.00U8 
-0.0055 



To estimate the sensitivity of this analysis without the effect of statistical correlations in the 
MC samples used, a toy MC study is conducted. In this study, angular distributions are generated 
for the signal and sideband regions to simulate the statistics available in the data and various 
MC samples used in the analysis. After subtracting background samples from the toy data, the 
angular distribution is corrected for efficiency and fitted using EqlH The statistical uncertainty 
on e obtained from the fit is 0.0063, which combined with the systematic uncertainties leads to an 
estimated uncertainty of (7^ =^00084. 

With the MC studies completed, the angular distribution in the data is obtained by subtracting 
estimated background events as described above. The remaining distribution is corrected for 
efficiency and fitted using Eq. [T]as shown in Figure |5l The fit has x^/dof = 15.4/18, and the 
fitted value of e in the data is (—5.5 ± 5.8(stat)l5;5(syst)) x 10~^, which is consistent with no SCC 



contribution to r 



U)1T 



Vt decays. For the upper limit on N', 



(non-vector current) ' (vector current) 



IK 



Bayesian approach |IT3| is used as negative values of e are non-physical. Using only the positive 
portion of the probability distribution for the value of the SCC contribution, etTue, the distribution 
for etrue is a bifurcated Gaussian with mean e = 5.5 x 10"'^ and errors a^ =l'^ooo84- The limits 
obtained from this method are 0.69% at 90% C.L. and 0.85% at 95% C.L. 
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Figure 5: The fitted cos 6^.,^ distribution for the data. The fitted curve is described in the text. 

5 SUMMARY 



A search for second-class currents in the decay r ^ tiJTr z/,- is conducted with the BaRar detector. 
No evidence for second-class currents is observed, and a 90% confidence level Bayesian upper 
limit for N'^Zn-^^aor current) / ^(vector current) ^^ ^^t at 0.69%. This limit is an Order of magnitude lower 
than the limit set by the CLEO collaboration ||3l. 
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